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Abstract. We study a stochastic predator-prey model on a square lattice, where 
each of the six species has two superior and two inferior partners. The invasion 
probabilities between species depend on the predator-prey pair and are supplemented 
by Gaussian noise. Conditions are identified that warrant the largest impact of noise on 
the evolutionary process, and the results of Monte Carlo simulations are qualitatively 
reproduced by a four-point cluster dynamical mean-field approximation. The observed 
noise-guided evolution is deeply routed in short-range spatial correlations, which 
is supported by simulations on other host lattice topologies. Our findings are 
conceptually related to the coherence resonance phenomenon in dynamical systems 
via the mechanism of threshold duality. We also show that the introduced concept 
of noise-guided evolution via the exploitation of threshold duality is not limited to 
predator-prey cyclical interactions, but may apply to models of evolutionary game 
theory as well, thus indicating its applicability in several different fields of research. 
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1. Introduction 

Cyclical interactions, the simplest non-trivial example being the children's rock-scissors- 
paper game, are despite their apparent simplicity fascinating examples of evolutionary 
processes [H E]. Examples from real life are the mating strategies of side-blotched 
lizards [3] , overgrowths by marine sessile organisms [1] , and competitions among different 
strains of bacteriocin-producing bacteria [5]. Moreover, studies of cyclical interactions 
may provide insights into the formation of defensive alliances |6l |7] and Darwinian 
selection |8], as well as structural complexity and pre-biotic evolution [9]. Thus, cyclical 
interactions are rightfully acquiring a central role in the study of evolutionary processes 
[To], not just in Lotka-Volterra like predator-prey systems, but also in evolutionary game 
theory [11] , where strategic complexity [121 US] often leads to a closed loop of dominance 
between participating strategies. 

An interesting concept that has only recently begun to supplement models of 
evolutionary processes is the addition of stochasticity at some level of interactions. While 
the origin of stochasticity can sometimes be related to the finiteness in population size 
[141 [T5] , uncertainties may also enter under the assumption of irrationality and errors 
in decision-making [161 [13 [IE]- Moreover, the approach is viable when stochasticity 
in payoffs originates from the heterogeneity of the system [19|, or if explicit payoff 
fluctuations are considered [20] . The latter examples are closely related to the concept we 
adopt presently, as stochasticity does not originate from finite populations but is due to 
uncertainties arising by predator-prey interactions. Irrespective of these particularities, 
however, stochasticity appears to posses the ability of having a profound impact on the 
evolutionary process [21] . Examples range from stochastic gain in population dynamics 
[22] and eradication of coexistence in the cyclic Lotka-Volterra model [23] , to cooperation 
promotion in the spatial prisoner's dilemma game [21] . 

Already long before making its debut by evolutionary models, stochasticity has 
been identified as a potentially crucial agonist in certain tj^es of dynamical systems. 
The so-called stochastic resonance, standing for the resonant noisy enhancement of the 
correlation between the system's response and a weak external stimulus, has been first 
reported for bistable systems [25], and latter on also for a broad variety of other physical 
as well as biological systems [26]. Importantly, noise can play an ordering role even in 
the absence of additional external signals, whereby the established term describing the 
phenomenon is coherence resonance [271 EHl EHl [301 EH [32] . Over the years a particular 
property of dynamical systems, termed excitability, has crystallized as being very 
beneficial when noise-induced phenomena are a desirable feature. Excitability, uniquely 
comprising elements of slow and fast dynamics thus enabling weak perturbations of 
the system to result in large-amplitude deviations before recovery, has fueled studies 
reporting various effects of noise on temporal and spatially extended dynamical systems 
for over a decade [331 [Ml [35] . It is the different noise dependencies of the slow and fast 
dynamics, constituting a so-called threshold duality, that are responsible for the majority 
of noise-induced phenomena observed in excitable systems [36]. In particular, while the 
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slow dynamics, representing the lower threshold, is very susceptible even to weak noisy 
perturbations, the fast dynamics, representing the upper threshold, is not. On the other 
hand however, the fast dynamics, resulting in large amplitude excitations, is sensitive to 
strong noise. Consequently, the temporal order of these excitations exhibits a maximum 
when we gradually increase the strength of noise [29]. 

Although some correlations between the coherence resonance in dynamical systems 
and noise-guided evolution have recently been established in terms of proximity to 
special bifurcation points [37j, the search for additional conceptual similarities linking 
the two avenues of research is still vibrant. This paper links these avenues by 
showing that a reticulate six-species predator-prey model with heterogeneous invasion 
probabilities offers adjustment possibilities via fine-tuning of a single parameter that 
takes the system from a noise frigid to a very susceptible state. In the latter state, 
the studied predator-prey model incorporates two thresholds that can be affected 
by noise, whereby one is small, conceptually similar to the slow dynamics, and the 
other large, conceptually similar to the fast dynamics of an excitable system. By 
introducing Gaussian noise, we show that this threshold duality can be readily exploited 
in a resonant manner, thus indicating conceptual similarities between the presently 
reported phenomenon of noise-guided evolution within cyclical interactions and that 
of coherence resonance reported previously for excitable dynamical systems. We 
emphasize, however, that the outlined similarities are indeed only conceptual, as the 
quantity presently of interest is not temporal order of the dynamics as by the observation 
of coherence resonance but the survival chance of different species within the habitat, 
which resonantly depends on the intensity of noise. Moreover, the threshold duality is 
not explicitly constituted by the slow and fast dynamics as by excitable systems, but as 
we will show below, is routed in the parameters defining the predator-prey interactions. 

The remainder of this paper is organized as follows. In the next Section we define 
the microscopic model and summarize the results for the noise-free case. We also 
describe the introduced Gaussian noise and comment on its impact on the elementary 
process. The details of Monte Carlo (MC) simulation are provided as well. In Sec. 3 we 
present the effects of Gaussian noise, and show that at appropriate system conditions 
the latter can have a large impact on the survival of different species, even reverting 
the evolutionary process in a resonant manner. In particular, while small and large 
intensities of noise fail to have a profound impact on the model, intermediate intensities 
induce a re-entrance effect of the seemingly defeated species, hence decisively guiding 
the evolutionary process. The non-monotonous effect of different noise intensities is 
supported by the application of the dynamical mean-field (DMF) approximation. To 
test the generality of our findings, we also perform MC simulations of the model on 
different host lattice topologies. At the end of Sec. 3, we show that the concept of 
threshold duality, linking the present results with those obtained previously for noise- 
driven excitable systems, is not limited to predator-prey cyclical interaction but can 
be applied to models of evolutionary game theory as well. Finally, Sec. 4 features a 
summary of main results and concludes by suggesting that the discovered mechanism 
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for noise-guided evolution is applicable under different circumstances, and could thus 
prove useful in various fields of research such as ecology or economy, where evolutionary 
processes are often subject to unpredictable factors. 

2. Predator-prey model 

The studied predator-prey model comprises six species that are initially uniformly 
distributed on the square lattice. The distribution function of species is given by a set of 
variables = 0, . . . , 5, where i runs over all L x L lattice sites. Predator-prey relations 
of nearest neighbors and the corresponding invasion probabilities (0 < a, (3,'y,6 < 1) are 
defined by the food web presented in Fig. [H More precisely, a determines the probability 
that an even labeled predator will invade an odd labeled prey, and vice versa for /3, while 
7 and 6 determine probabilities of invasions within even and odd species, respectively. 
As exemplified in [3HI EH] for lattice Lotka-Volterra type models, the studied six-species 
model can be defined by a concise reaction scheme. In particular, for species s = 2k 



s + j ^ 2s (1) 

s + I ^ 2s (2) 

and for species s = 2k + 1 

s + j A 2s (3) 

s + I ^ 2s (4) 



where j = {s + 1) mod 6, / = (s -|- 2) mod 6, and A; = 0, 1, 2. In case of homogeneous 
invasion probabilities {a = /3 = j = 6=l) the system has two equivalent three-species 
states [denoted by (0-1-2-1-4) and (1-1-3-1-5)] exhibiting a self-organizing pattern 
maintained by cyclic invasions. Departing from a randomly distributed state where 
every lattice site is occupied by one of the six species, a spatial organization of growing 
domains begins. The domains are formed either by the members of the red alliance 
(0 -|- 2 -|- 4) or the members of the blue alliance (1-1-3-1-5) (see Fig. [1] for the color 
assignment). In accordance with the coarsening process on the spatial grid, one of 
the two domains is eliminated with equal probability, and the system ends up in the 
absorbing state consisting of members of either the red or the blue alliance only. In 
order to characterize the stationary state more precisely, we may introduce the order 
parameter m = pi + Ps + P5 — Po — p2 — Pa, whereby ps (s = 0, . . . , 5) denotes the fraction 
of species s on the spatial grid. Here m = 1 corresponds to the exclusive presence of 
the blue alliance (1-1-3-1-5), while m = — 1 indicates the absolute authority of the red 
alliance (0-1-2-1-4). The two domains are more precisely called "defensive alliances" 
because their members protect each other cyclically against the external invaders P, [7j . 
For example, if species invades the allies (l-f 3 -f 5), in particular by attacking species 
1, its intention is immediately disabled by the species 5 that is superior to both and 1. 
Thus, the intruder is quickly abolished from the (1-1-3-1-5) domain by the very same 
species 5 that dominates species 1 within the alliance. This reasoning applies for all 
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Figure 1. Food web of the studied predator-prey model. Arrows point from predators 
towards prey with heterogeneous invasion probabihties specified along the edges. Inner 
loops are colored for easier terminology. 

other possible attempts of non-allied species to invade a defensive alliance. Noteworthy, 
in this mechanism the proper spatio-temporal distribution of species plays a crucial 
role, meaning that the classical mean-field approximation, assuming a well-mixed state, 
cannot reproduce this feature. An interesting possibility is now to study effects of 
inhomogeneous invasion probabilities, which might affect the evolutionary process in 
a non-trivial way [IQ], essentially forming the basis for noise-guided evolution to be 
reported below. 

Since effects of inhomogeneous invasion probabilities on the evolution of defensive 
alliances in the studied cyclically dominated model have already been presented in [H], 
we here just summarize the main findings that are essential for the present study, and 
refer the reader to the previous work for more details. First, it is important to note 
that heterogeneities in the invasion probabilities are introduced in an alliance-specific 
way. In particular, we can study what happens when one defensive alliance is more 
aggressive towards the other {a ^ p), or when the internal mechanism fails to assure 
flawless protection against the invaders (7 7^ 6). In order to address these two issues 
systematically, we introduce two control parameters that, due to symmetries in the 
food web, uniquely determine the stationary state of the system characterized by the 
order parameter m. Namely, let G = (3 — a and H = 'j — 6 where H,G & [—1, 1]. 
Since G < {/S < a) and H > > 6) clearly favor the survival of the red alliance 
(0 + 2 + 4) (and vice versa for the blue alliance), an interesting competition obviously 
emerges only if G > and if > (or equivalently if G < and H < 0). Due to the 
symmetry of the problem we restrain or study to the parameter space spanning over 
H,G & [0, 1]. As soon as G rises above zero (keeping H = 0) allies (1 + 3 + 5) are 
favored as their members invade non-allied species more successfully {a < P), thus MG 
simulations yielding m = 1 in the stationary state. However, the advantage of the blue 
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alliance (1 + 3 + 5), given by G > 0, may be compensated by choosing sufficiently large 
values of H. Especially, if H rises above zero the internal invasions within the alliance 
(1 + 3 + 5) slow down in comparison to (0 + 2 + 4), thus decreasing the effectiveness of the 
protection shield of the blue alliance and in turn nullifying its advantage given by G > 0. 
Surprisingly the results of MC simulations reveal a non-monotonous phase diagram in 
dependence on H, as shown in Fig. [2l In particular, the advantage of aUies (1 + 3 + 5) 
again increases if H approaches 1 (the internal invasion probability S vanishes). Note 
that as (5 — allies (1 + 3 + 5) essentially stop to invade each other within the alliance. 
Via a stability analysis of the interface separating members of the two defensive alliances 
we have shown [41j that the penetration of the red alliance into the unfavorable G > 
region is a consequence of an interface driven effect. Presently, we will show that the 
spatiality of the system is also a relevant feature when the system is driven from a noise 
frigid to a noise susceptible state by varying the parameter H. 

Turning to the effects of noise, we should point it out that the model involves 
stochasticity already in its present form. In particular, the invasions of predators 
to neighboring prey sites are not deterministic but characterized by a probability. 
The latter source of stochasticity is routed in the spatiality of the model, dictating 
a random selection of one neighbor of each species in order to carry out an elementary 
invasion process. Stochastic effects also cannot be neglected if the system size is finite. 
However, this is a plausible assumption for populations studied within the context of 
evolutionary game theory [191 1121 US] . In these systems the strategy adoption is also 
ruled by a probability function that depends on the payoff differences between the 
two competing strategies [H]. Motivated by bounded rationahty, often an additional 
source of stochasticity is introduced by the so-called smooth imitation rule that allows an 
inferior strategy to replace a more successful one. The application of the Fermi function, 
adopted from statistical physics, allows to measure the intensity of selection via a single 
parameter K, termed accordingly as the temperature of selection. In the K 
limit the stochastic effects can be neglected, while for X oo limit the stochasticity 
is maximal. In our predator-prey model, lacking individual payoffs, the stochasticity 
resulting from finite K can be modeled by introducing Gaussian noise additively to 
all invasion probabilities so that the modified invasion probability is p = x + ^ where 
X G {a,(3,-y,6) and ^ is a random variable being normally distributed with zero mean 
and a standard deviation. It is worth noticing that the introduced Gaussian noise is 
in fact multiplicatively coupled with the dynamics of ps, as can be inferred from the 
classical mean-field equations or the DMF approximation. The Gaussian noise satisfies 
the correlation function 

where indexes m and n mark the central location of two predator-prey pairs on the 
lattice, while h and k denote two consecutive pair interactions. Thus, the stochastic 
disturbances are delta correlated (uncorrelated) in time as well as among predator-prey 
pairs. Importantly, if an invasion probability becomes negative {p < 0) due to noise we 
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Figure 2. Phase diagrams in dependence on tlie two main system parameters H and 
G for various a. Lines denote the border separating the two pure phases that are 
characterized either by the exclusive dominance of the red (0 + 2 + 4) (m = — 1) or the 
blue (1 + 3 + 5) (m = 1) alliance. 

allow backwards invasions (a prey can occupy a predator's site), while potential values 
of invasion probabilities above 1 {\p\ > 1) are simply treated as sure events. We thus 
allow noise-induced reversals of the direction of dominance depicted in Fig. [H 

Each elementary step of the MC simulation involves two actions. First, two nearest 
neighbors are chosen at random, and second, if the two neighboring species form a 
predator-prey pair (species directly connected by an arrow in Fig. [1]) the prey (or 
predator depending on ^) is killed with the probability p = x + ^, where x is the 
invasion probability specified along the connecting arrow and ^ is the Gaussian noise 
with properties as described above. On the other hand, if the two randomly chosen 
species form a neutral pair (species not directly connected by an arrow in Fig. [1]), or 
if both are identical, the second part of the elementary MC step dictates no action. 
In accordance with the random sequential update, each individual is selected once on 
average during a particular MC step. 

3. Results 

First, we study how the phase separation line in the H — G parameter space varies in 
dependence on a. Due to symmetries in the food web, and to ensure transparency, we 
perform all below calculations using /3 = 1.0 and 7 = 1.0, thus designating a and 6 as 
the only parameters able to change H and G. Figure [2] shows the results, from which it 
can be inferred that there exist an intermediate a for which allies (0 + 2 + 4) receive the 
biggest boost, or in other words, are able to compensate for the largest values of G by a 
given H. Clearly, there exist a resonant dependence on a as larger levels of stochasticity 
fail to have the same effect because the phase separation hne gradually converges to the 
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Figure 3. Absolute rise of the phase separation line AG in dependence on a for 
different values of H (pink= 0.1, blue= 0.2, green= 0.5, red= 0.9). Note that positive 
values of AG manifest a facilitative effect of noise on the survival of the red alliance 
(0 + 2 + 4). The inset shows results from the four-point DMF approximation for the 
same values of H as in the main figure, in particular depicting critical values of G 
separating the two pure phases. 



zero-noise limit as a increases. 

To study the outlined resonant dependence, we plot AG as a function of a, 
whereby AG measures the rise of the phase separation line in the vertical direction 
with respect to its position at cr = at any given H. Results for different values of 
H are presented in Fig. [3l Two important conclusions can be made. First, added 
Gaussian disturbances affect the evolutionary process resonantly in dependence on a, 
and second, the phenomenon becomes increasingly pronounced as H ^ 1 (6^0). We 
should also point out that by values of H close to noise has absolutely no impact on 
the evolutionary process. Thus, we demonstrate that as the parameter H is varied from 
to 1 the system goes from a noise frigid to a noise susceptible state, and moreover, 
noise-guided evolution sets in if a is appropriately adjusted. 

Our numerical findings can be supported by an alternative approach entailing 
the application of /c-site cluster DMF approximation. Since the classical mean-field 
approximation assumes well-mixed populations, it is easy to understand that it cannot 
reproduce above-described results. As we have argued earlier, the spatiality plays a 
fundamental role by the explanation of the penetration of the red alliance into the 
unfavorable G > area. The shortage of the classical mean-field approximation 
may be eliminated by applying the DMF approximation technique that proved to be 
very appropriate for obtaining qualitatively correct phase diagrams for several non- 
equilibrium systems. For a detailed description of this method we refer the reader to 
earlier papers HH] . The DMF approximation is a dynamical version of the cluster- 
variational method, and it involves finding a hierarchy of evolution equations for the 
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probability distributions of configurations within a cluster of k sites. If correlations 
in larger clusters are neglected then the level of approximation can be characterized 
by the value of k. Similarly as by the MC simulations, in the DMF approach the 
introduction of Gaussian noise can, besides altering the strength of invasions, also revert 
the direction of dominance. Therefore, aside from handling the elementary processes 
displayed in Fig. [1], we also have to consider the possibihty of reverse processes. With 
this extension, the four-point level of approximation can describe correctly not only the 
impact of the dynamical benefit of the red alliance, but also the impact of Gaussian noise. 
Indeed, displayed as an inset in Fig. [3l the re-entrance phenomenon driven by increasing 
values of a is well reproduced, and moreover, exhibits the same dependence on H as 
revealed by MC simulations. In particular, values of H close to zero yield a virtually 
noise-resistant dynamics across the whole span of a, whereas noise-guided evolution 
sets in as — > 1 for intermediate levels of stochasticity. Nonetheless, a slight difference 
between the results of MC simulations and the DMF approximation can be observed. 
According to the latter, large values of a can eliminate the dynamical benefit of allies 
(0 + 2 + 4), thus resulting in the solution predicted by the classical mean-field theory. 
An even higher level of approximation might still eliminate this minor inconsistency 
with the MC simulations, but unfortunately the resulting number of variables is too 
large to be manageable by the computer resources currently available to us. In sum, 
the failure of classical noise-driven mean-field equations to reproduce the results of MC 
simulations and the necessity for the four-point cluster DMF approximation confirm 
that the reported noise-guided evolution within cyclical interactions is indeed deeply 
routed in the short-range correlations of the spatial distribution, as conjectured already 
from the outlay of phase diagrams depicted in Fig. [2] and findings reported previously 
in [H]. 

In order to strengthen our arguments, we also examine the dynamics of the noise- 
driven model via MC simulations on two additional types of regular graphs, namely on 
the honeycomb and the triangular lattice. Our choice is motivated by the fact that the 
two types of host lattices have different coordination number z. While the honeycomb 
lattice {z = 3) allows manifesting the importance of spatial correlations, the triangular 
lattice {z = 6) may move the system closer towards mean-field conditions. Therefore 
we argue that the noise driven re-entrance phenomenon will be more (less) pronounced 
when the honeycomb (triangular) host topology is applied. Indeed, results presented in 
Fig. m fully support this expectation as the observed phenomenon is substantially better 
expressed for z = 3 than for z = 6, and thus additionally strengthen the fact that the 
reported noise-guided evolution within cyclical interactions is routed in the short-range 
correlations of the spatial distribution. 

Next, it is of interest to draw some analogies between the reported noise-guided 
evolution within cyclical interactions and the phenomenon of coherence resonance 
reported previously in excitable dynamical systems [291 EH]. Although we presently 
consider the survival of different alliances as being resonantly dependent on a and not 
the temporal order of their evolution, and are thus strictly speaking unable to write 
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Figure 4. Absolute rise of the phase separation line AG in dependence on a for 
the honeycomb (green) and the triangular (red) lattice, obtained by _ff = 0.99. The 
interpretation of AG is the same as by Fig. [3l Importantly, the attainable effect of 
noise decreases as the coordination number z increases. 



about a classical coherence resonance phenomenon, we still argue a conceptual link 
can be established by addressing the presently reported phenomenon as an evolutionary 
coherence resonance [3^. Noteworthy, related conceptual differences emerge also in |22j . 
where the authors employs similar arguments to address the problem. Aside from this 
difference, however, we will show that the threshold duality, constituted by the noise 
susceptible slow dynamics and the noise robust fast dynamics of excitable systems, can 
be identified also in the present model, and is in fact crucial for the observation of 
noise-guided evolution. 

The threshold duality can be made visible by studying the properties of invasion 
probabilities along the phase separation line of the deterministic model shown in Fig. [2j 
We start with close to zero values of H where, taking into account results presented in 
Fig. [31 noise has no impact on the evolutionary process irrespective of a. There, the 
switch between the dominance of allies (1 + 3 + 5) and (0 + 2 + 4) occurs at = 0.1 and 
G = 0.065, or equivalently, at 6 = 0.9 and a = 0.935 (note that f3 = 1.0 and 7 = 1.0 
are held fixed throughout this work). Importantly, the differences 6th = 1 — and 
ath = 1 — a define the two thresholds Gaussian noise can influence in order to guide the 
evolutionary process. In particular, if one of the two invasion probabilities a or 6 would 
become equal to 1, the corresponding defensive alliance would benefit and eventually 
defeat the other, as can be inferred from the read line depicted in Fig. [2] (consider for 
example the case when 6 = 1.0 while a = 0.935; clearly the blue alliance wins making 
m = 1 in the stationary state). However, since at H = 0.1 both thresholds have virtually 
identical magnitude {6th = 0.1 and ath = 0.065), and are exposed to the same intensity 
of temporally and spatially uncorrelated Gaussian noise, the long-term probability of 
exceeding one threshold while leaving the other untouched is practically zero. Thus, 
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noise cannot influence the evolution of the two defensive alliances irrespective of a, as 
displayed in Fig. [31 However, the situation changes substantially as H increases. By 
H = 0.5 the phase separation line is situated at G = 0.175 in the noise-free case. The 
two relevant thresholds then equal 6th = 0.5 and ath = 0.175. As above, if Gaussian 
noise would exceed only ath the red alliance would win (m = —1), and via similar 
reasoning, if only dth would be exceeded the blue alliance would win (m = 1). Contrary 
to the H = 0.1 case here the two relevant thresholds differ substantially. Therefore it is 
easy to see that by an appropriate value of a the threshold ath = 0.175 will be exceeded 
statistically far more often than 6th = 0.5, which ultimately results in noise-guided 
evolution, specifically giving the seemingly defeated alliance (0 + 2 + 4) the winning 
edge over the allies (1 + 3 + 5), as indicated by the positive values of AG in Fig. [31 
However, by larger values of a the upper threshold 6th = 0.5 will be exceeded often as 
well, thus again nullifying the dynamical benefit coming from the faster inner cycle for 
the red alliance, which explains the depicted resonant dependence of AG. Note that 
although larger a increase the probability of crossing the lower threshold {ath = 0.175) 
as well, the maximal value of each invasion probability equaling 1 (above that each 
event is simply treated as a sure event) prohibits a noticeable effect of this fact, and 
thus results in the decline of AG as a increases. We argue that the difference in the two 
thresholds (lower threshold given by ath = 0.175 and the upper one given by 6th = 0.5), 
constituting a threshold duality in the presently studied model of cychcal interactions, 
is conceptually similar to the threshold duality observed in excitable dynamical systems, 
albeit the latter is routed in the fast and slow dynamics of the system [29l [36] , while the 
former originates directly from the parameters values defining the evolutionary process. 

Finally, it is straightforward to extend the above argument also to the case when 
H = 0.99, whereby the noise-free phase separation line is situated at G = 0.09, and hence 
the two relevant thresholds, equaling 6th = 0.99 and ath = 0.09, clearly differ by an order 
of magnitude thus constituting the vital threshold duality needed for the observation 
of noise-guided evolution. Importantly, an additional threshold emerges related to 
the possibility of noise-induced negative values of 6. Specifically, the new threshold 
simply equals to 6 = 0.01, and is apparently of similar magnitude as ath = 0.09. Since 
negative values of 6 flip the circle of dominance within the blue alliance (1 + 3 + 5), thus 
temporarily disrupting its protective shield against the invaders, intermediate values 
of a warrant a two-fold advantage to the red alliance, ultimately resulting in the best- 
expressed resonance curve in Fig. [3l The origin of the two- fold advantage of the red allies 
can be studied more precisely by measuring the success rate of invasions along the two 
inner circles (depicted with blue and red in Fig. [1]) of the food web. To be more specific, 
we calculate the normalized difference of clockwise and anti-clockwise invasions w for 
the red and blue alliance in dependence on a. Note that it; = 1 if all invasions occur in 
the anti-clockwise direction (as depicted by arrows in Fig. [1]) and zero if clockwise and 
anti-clockwise invasions occur equally often. In other words, the value of w characterizes 
the strength of the defensive mechanism within an alliance. Results in Fig. [5l show that 
the defensive shield of the blue alliance (1 + 3 + 5) suffers virtually immediately after 
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Figure 5. Normalized difFerence of clockwise and anti-clockwise invasions w along the 
two inner circles (depicted with blue and red lines in accordance with the coloring in 
Fig. [T]) of the food web in dependence on a. Panel (a) show results for H = 0.99 and 
G = 0.07, while panel (b) features results for _ff = 0.1 and G = 0.09. 

noise is introduced, whereas the red alhance is able to uphold a perfectly functioning 
protection until a ^ 0.6. Naturally, as a increases further the defensive capacity of 
the red alliance also weakens. However, the window resulting from the delayed faint 
of the defensive mechanism, in conjunction with the previously described mechanism 
warranted by the threshold duality via 5th and ath, enables a resonant strength of the 
referred (0 + 2 + 4) domain. Results displayed in Fig. [3] support this argument since 
a ~ 0.6 limits the peak of the resonance curve obtained for H = 0.99, thus signaling 
the beginning of end of the noisy support for the red allies. Also in accordance with the 
results presented in Fig. [3l for H = 0.1 the dependence of w on a is very similar for both 
defensive alliances (red and blue), thus clearly signaling the absence of the mechanism 
warranting noise-guided evolution in the studied model. The simultaneous decrease 
of the defensive strength results in a noise-frigid state since all invasion probabilities 
are of similar magnitude. The above-described simple arguments provide a workable 
mechanism for noise-guided evolution, and at least conceptually unify the phenomenon 
with the coherence resonance in excitable dynamical systems. 

To end this section, we would like to note that the above explanation warranting 
the resonant dependence of AG on a, intimately relying on the existence of threshold 
duality in evolutionary processes, is valid for other systems as well. To demonstrate this, 
we briefly consider a spatial prisoner's dilemma game [IHl HTJ HH] , where two cooperators 
receive the reward R = 1, two defectors receive the punishment P = 0, whilst a 
cooperator and defector receive the suckers payoff S = —r and the temptation T = 1 + r, 
respectively, thus satisfying the prisoner's dilemma payoff ranking T>R>P>Sif 
r > 0. It is easy to identify a threshold duality in the system if only r ^ 1. If 
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the payoffs become stochastic as a result of additive Gaussian noise JIH], inequahties 
T > R and P > S differing by r will be violated statistically much more frequently 
by a given a than the R > P inequality differing by 1. It is easy to understand that 
intermediate levels of stochasticity promote cooperation because the violation oiT > R 
and P > S favors the cooperative strategy since it potentially nullifies the advantage r 
defectors have over cooperators (two cooperators might end up receiving a larger payoff 
each than a defector facing a cooperator, and a cooperator facing a defector might be 
better off than two defectors). On the other hand, this facilitative effect is limited by 
violations of R> P, belonging to the so-called social dilemma [SOlET], because then two 
defectors might be better off than two cooperators, which again gives the winning edge 
to the defecting strategy, and hence results in a resonant dependence of cooperation 
fitness on a |20j. Although being fairly simple, the described explanation outlines a 
general mechanism of cooperation promotion in the spatial prisoner's dilemma game, 
thus giving another example of noise-guided evolution relying on the threshold duality 
in evolutionary processes. 

4. Summary 

In sum, we provide conclusive evidences that the threshold duality is crucial for the 
observation of noise-induced resonances in evolutionary processes. Particularly, we 
show that a reticulate six-species predator-prey model with heterogeneous invasion 
probabilities possesses a single parameter that is able to guide the system from a 
noise frigid to a noise susceptible state, thereby relying exclusively on the emergence 
of threshold duality. The phenomenon revealed by MC simulations can be verified 
by a four-point cluster DMF approximation, thus providing ample evidences that the 
reported noise-driven evolution within cyclical interactions is heavily routed in the short- 
range correlation of the spatial distribution. This conjecture is additionally supported 
by studying the evolution of species on different regular graphs, and confirmed by the 
fact that the phenomenon of noise-guided evolution weakens as the coordination number 
of the host lattice increases. While the observations in the studied model are related 
to the system-specific non-monotonous phase separation line between the two defensive 
alliances, we note that the general concept is valid also in the framework of evolutionary 
game theory, thus suggesting it may be widely applicable in various fields of research 
ranging from economy to ecology. 
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